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Abstract 

A companion article analyzed very weakly first-order phase transitions in 
the cubic anisotropy model using e expansion techniques. We extend that 
analysis to a calculation of the relative discontinuity of specific heat across 
the transition. 
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I. INTRODUCTION 



In this work, we use e expansion methods to compute the universal ratio C + / C_ of specific 
heats for arbitrarily weak first-order phase transitions in the cubic anisotropy model. C + and 
C- are the specific heats of the disordered and ordered phases at the transition temperature. 
This work is a direct follow up to the computations of other universal ratios in ref. and 
we shall eschew any discussion of motivation or review of method and notation; instead, we 
jump directly into the calculation. The reader should read ref. |I] first. 

In the next section, we review the leading-order calculation of the ratio, which was first 



performed by Rudnick in ref. ||. Our result differs by a factor of 4. In section |T|, we 



proceed to next-to-leading order in e. Our final results are displayed in section IV 



II. LEADING-ORDER ANALYSIS OF C+/C 

Recall from the introduction of ref. ffl] that, in the three-dimensional theory, the square 
m 2 of the scalar mass plays the roll of the reduced temperature near the transition and 
the effective potential V represents the free energy of the system. The specific heat can be 
extracted from the effective potential as 

The proportionality constant will not matter to our determination of the ratio C + /C_.0 

Because of this relationship, it is useful to begin by summarizing the leading-order form 
of the effective potential discussed in ref. IJ. 



^dn particular, the reduced temperature is proportional to m 2 (/io) for some fixed renormalization scale 
but we will instead apply ([!]) with m 2 (fi) where ufi is roughly the order of the correlation length and varies 
as we approach the transition. [Specifically, we choose /i so that u(fi) = — v(fi).] However, m 2 (^t) is related 
to m 2 (/^o) by multiplicative renormalization which, even though /i dependent, cancels in the ratio. 
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A. Summary of leading-order potential 



At one-loop order, the effective potential along an edge 
the tree-level instability line u— — v, is [eq. (3.9) of ref. ||T|]: 
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where M 2 = j;Afu(j) 2 , Cn = Mn — 1), Ci = — §Cn, and the normalization A/" is 



0( 



(2) 



A/"= (47r) d/2 r(f-i) 



(3) 



in d = 4— e spatial dimensions. 

A first-order transition occurs as m 2 is varied. At the transition, in the asymmetric 
phase, m 2 /M 2 ~ 0(e). In the asymmetric phase, the potential then reduces to 

A//i e (Vo + Vi) = A + 3u~ l m 2 M 2 + M 4 C n In f — J + C 10 + 0(e 2 V% ym . (4) 

The transition occurs when 

m 2 = m \ [1 + 0(e)] , 
and the asymmetric minimum is at 

M 2 = M 2 [l + 0(e)] , 



(5) 



(6) 



where 



2 CuUfi 2 ( C 10 \ n-l 2 1/2 
M 2 ^ 2 exp ("1-^)=A V2 . 



(7) 
(8) 



We will also need the variation of the asymmetric minimum as one varies m 2 slightly away 
from the transition: 

d(M 2 ) 3 



d(m 2 



i 2 = m 2(l+0(e)) C U U 
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+ 0(e c 



(9) 



Finally, we will need the renormalization group flow of the various couplings, previously 
presented in ref. [[[]. At leading order, 

^m 2 = (^l(u,v) + 0(u 2 ,v 2 ))m\ (10) 



fid^u = -eu + (3^\u,v) + 0(u 3 ,v 3 ) 
/i<V' = -eu + Pl 1] (u,v) + 0(u 3 ,v 3 ) . 



where 



( i) _ n + 2 



u + v, 



/?(!) = v{Au + 3v) . 
It is also convenient to introduce / = u/v and 

^f = pf\uj) + 0(u 2 ,v 2 ) 
n - 4 



(11a) 
(lib) 



(12) 
(13) 
(14) 



(15) 
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B. Specific heats 

The most non-trivial task in computing the specific heat ratio at the transition will be 
to handle the renormalization group flow of the constant term A in the potential. Before 
investigating this, we first write down the expression for C + and C_ given the potential. 
From (0), 
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(17) 



where the second derivative of A is understood to be evaluated at the transition m? = 
m\(l + 0(e)). Using (|j) for the asymmetric phase, and remembering that 8mV(M) = at 
the asymmetric minimum, 



d ^ 2 



m 2 = m 2(l+0(£)) 



d(m 2 



d \ 2 . 3d(M 
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A + - 



d(m 2 ) J u d(m 
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+ 0(l/e) 



-\d {m 2 )) A -^ + 0(1/6) ' (18) 
where M has been evaluated at the transition using (P) and ©. C + will be shown to be 
0(1/ 1) while C_ is 0(l/e 2 ). 

C. The running of A 

Now we turn to the contribution to the specific heat from the term A. The renormaliza- 
tion group equation for A at one-loop order is 

Ti 

nd fl A = eA + -m\l + 0(u,v)). (19) 
The solution, in terms of m 2 (/z), is 

m = r^t^Xim^rii+oM), ( 20 ) 

y/i / ^ Jfj, /i \/i / 

where fio is some initial scale. The running (|T0|) of m 2 yields 

m 2 (/i) = E^\fi, /i ) m 2 (/i ) (1 + 0(u, «)) , (21) 

where 

i^Wo)=exp(f (22) 

and the integral is to be understood as evaluated along the leading-order solution for the 
coupling constant trajectory. The contribution of A to the specific heat is then 
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d 



d(m 2 



A = n 



^(^] e [E^^,^}\l + 0(u,v)). 



(23) 



can be easily evaluated by changing variables from \j! to / using ( p"5"D and noting that 
tftf/P? (OS) depends only on /: 

(raA-l)/2 
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where 
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(25) 



The remaining integral in (|23|) can be performed similarly if we make use of the following 
relation ((AlOa) of ref. 0) for the solution to the leading-order RG flow equations (|TT1): 



wo ffo 



(26) 



This gives 
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nA 



df 



/; v/+a 



(l + 0(u,v)) 



u(f + A) J/o 



So the result is 



d 



d(m 2 



A 



df(l + 0(u,v)) 



n\ [/(/i) - /(AH>)] 



(27) 



[1+0(6)]. 



(28) 



u{n) (/(/,) + A) 

The only other elements we need are the values of f(fio), /(/•*), and w(/i) for the desired 
trajectory. As discussed in ref. HJ, we can obtain the universal ratios of interest by studying 
the trajectory that flows away from the cubic fixed point at /(/io) — ~A to the line of 
classical instability u— — v at /(//) = —1 and u = it*, where, at leading order, 

3(n 2 + 5n + 3) 



n(n + 2)(n + 



-e . 



(29) 



For this trajectory, 



This is 0(1/ e) and dominates C + ([TTD , but it does not contribute at leading order to C_ 
(|T8|). Putting it all together, 

C+ nXC u 
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e + 0(e 2 ) for n = 2. 



(31) 



This result is 4 times smaller than the result originally quoted by Rudnick [fj]. 

III. NEXT-TO-LEADING-ORDER ANALYSIS OF C+/C 

Because the contribution of A dominated C + , our formula ( |TT| ) for C + is adequate at 
next-to-leading order provided we extend our analysis of A to next-to-leading order. Before 
doing so, let us first consider C_. 



A. The asymmetric phase: C_ 

The two-loop effective potential near the asymmetric phase is given in ref. 
eqs. (4.1), (4.2), and (4.4)): 

A///(V + V x + V 2 ) = A + Zu-^M 2 + M 4 



where 
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As discussed in ref. U, we shall not need to know C e i and C 20 . The transition takes place 
at 



2 2 

m = m 2 



l + 0(e 2 



with the asymmetric minimum at 
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At the transition, 
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The result for C_ is 
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d(m 2 ) J u 2 Cu 

Since the contribution of A is sub-leading, the leading-order result ( |3"0"D for the contribution 
is adequate here. 



B. The NLO running of A 



Now we are left with calculating d 2 A/d(m 2 ) 2 to next-to-leading order. The two-loop 
renormalization group equation is 
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?? 



fid.A = eA + ^(m 2 ) 2 (l + 0(u 2 , v 2 )) 



(42) 



which has the same form as the one-loop equation. It's solution then also has the same form, 



'sr)) A=»jr^(^)'w,rtp(i + o(«',o) 



(43) 



where 



E(ii, no) = exp ( f ^-p m 2 J 



(44) 



but now both the /^-functions and renormalization group trajectories should be evaluated 
at two loops. The two- loop /^-functions are given in Sec. VIA of ref. M: 



p = pm+pw + o(u 4 ,v 4 ), 



(45) 



where 
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To evaluate the integrals, we again change variables from \x to /, and we shall treat 
the two-loop effects on /3-functions and trajectories perturbatively. Following sec. VI D of 
ref. [Q, it is helpful to make the e dependence explicit by rewriting (it, v) = e(u, v), and the 
expansion of the trajectory gives 



u(f)=u^(f) + e5(f) + 0(e 2 ), 

where u^(f) = fR{f, c) is the one-loop result described in sec. IV C of ref. 
for 8(f) is given by eq. (6.23) of ref. [|I|]. 
To change variables from ji to /, we use 



(49) 



The solution 



df 



P 
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1 - 



«(/),/ 

(5{f)d a pf + pf) + 0{<? 



(50) 
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The subscript il w^(f), /" at the end of this equation means that the /3-functions in the 
expression are to be evaluated with u — > u^(f) and v — > u^(f)/ f. Expanding the definition 
((HI) in e then gives 



E(f, f ) = f ) [1 + e6E(J, /„) + 0(e 2 



(51) 



where E^ is the leading-order form of ( f2"4"| ) and 



Pm2 



(1)^ 



(52) 



aW (/')>/' 

To do the final integral of (|43"D, we need an expansion of (/x//x') £ . This can be obtained 
by writing 



y =exp V V ~W , 



(53) 



and converting to / with the expansion (p0| ). Putting all the expansions together yields 

2 
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where 



X 1 (/') = --| I j(W^/3} 1) +/3f 
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This integration does not seem to have a simple form for general n. For n = 2, we are able 
to obtain a simple result for the trajectory flowing away from the cubic fixed point to the 
classical instability line u = —v: 



d 



d(m 2 



A = U + e 

u, 



49 3 3 

In - 

40 5 2 



+ 0(e 2 )} 



(57) 



11 



Putting everything together, and using the next-to-leading order n=2 result 



51 
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(58) 



for u* from eq. (6.28) of ref. ffl, our final result is then 



C + 17 r 17 , 
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(59) 



for n=2. 



IV. DISCUSSION 



Evaluated numerically, the final result ( 59f) for the ratio is 

a 



— - = 0.0531e 
c _ 



1 + e(-0.2125 In e + 1.3993) + 0(e 2 ) . (60) 



This ratio is compared against Monte Carlo simulations in ref. [|J. The 140% correction 
at next-to- leading order for e—1 suggests that the e expansion will be at best marginally 
successful for this quantity. 
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